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Abstract 

Radio variability on timescales from a few hours to several days in extra- 
galactic flat-spectrum radio sources is generally classified as intra-day variability 
(IDV). The origin of this short term variability is still controversial and both 
extrinsic and intrinsic mechanisms must be considered and may both contribute 
to the observed variations. The measured linear and circular polarization of IDV 
sources constrains the low energy end of the electron population. Any popula¬ 
tion of cold electrons within sources at or above the equipartition temperature of 
10 11 K depolarizes the emission and can be ruled out. Intrinsic shock models are 
shown to either violate the large fraction of sources displaying IDV or they do 
not relax the light travel time argument for intrinsic variations. From structure 
function analysis, we further conclude that interstellar scintillation also leads to 
tight size estimates unless a very local cloud in the ISM is responsible for IDV. 
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1 Introduction 

Intraday Variability (IDV) of flat-spectrnm compact Quasar cores and BL Lacs at 
cm-wavelength has been discovered in 1985 (Heeschen et al., 1987) and is a common 
phenomenon (~ 25%) (Quirrenbach et ah, 1992) among these sources. On rare occa¬ 
sions in 0716+714 and 0954+658, correlations with optical variations have been found 
and the evidence has been reviewed several times (e.g. Wagner S.J., Witzel A., 1995; 
Krichbaum, this volume). Radio-optical correlations, if real, suggest either fast in¬ 
trinsic variations or gravitational lensing. From the light travel time argument the 
apparent brightness temperatures A are in the range of 10 17 -10 21 K and far in excess of 
the intrinsic inverse Compton (IC) limit of 10 12 K. The observed superluminal motions 
of jet components in many of these sources imply T = 5-10 and allow for Doppler 
boosting factors V = (T(l — f3 cos O))^ 1 an d time shortening, which are insufficient for 
reducing the apparent A down to the IC limit. The required Doppler factors V are 
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in the range of 60 up to 1000. Furthermore jets with a surface brightness at the IC 
limit are radiatively inefficient (Begelman, Rees & Sikora, 1994) and carry most of 
their energy as bulk motion. This is not supported by the observed power of radio 
lobes of these sources and raises the energy requirement for the central engine to an 
uncomfortable level. It is therefore argued (Readhead, 1994; Begelman, this volume) 
that incoherent synchrotron sources in jets should radiate at the equipartition temper¬ 
ature Te ~ 10 11 K. This enforces Doppler factors which are 2-3 times larger than for 
the IC limit. Furthermore synchrotron sources at the IC limit with bulk motions of 
T > 100 are dominated by inverse Compton scattering of either AGN photons (Begel¬ 
man, Rees & Sikora, 1994) or CMB photons at redshifts z ~ 1 and not by the SSC 
process. The cooling is catastrophic, independent of the brightness temperature, and 
this explanation must therefore be discarded. 

Other more tempting suggestions are the propagation of relativistic thin shocks in 
the jets (Qian et ah, 1991), so that the observed variability timescale is not a measure of 
the source size, and scintillation induced by the interstellar medium (e.g. Rickett 1990) 
of otherwise non-variable sources. We will explore both explanations in the following 
sections. 


2 Cool Particle Depolarization 


The emission of incoherent synchrotron sources with high brightness temperatures is 
dominated by radiation from the r = 1 surface. For intrinsic TVs above the equiparti¬ 
tion temperature, the energy of the source is dominated by particles (e~, e + ,p) with a 
strong dependence on T b : 

U e /U B (x T 7+2a , 

where a is the optically thin spectral index of the emission. Above the IC limit any con¬ 
tribution of cold electrons in an e“/p plasma will depolarize the synchrotron emission. 
Emission from the r = 1 surfaces of flat spectrum radio cores comes predominantly from 
electrons with constant 7 ra d-factor, independent of frequency. We consider power-law 
distributions of electrons N( 7 ) oc 7 -p above a lower cut-off 7 min . The degree of linear 
polarization is likely to be reduced by tangled magnetic fields in a turbulent plasma 
with a typical wavenumber /c 0 in sources of size R (e.g. jet radius). For independent 
variations of magnetic held orientation in adjacent cells of size /cj 1 the fractional linear 
polarization is 717 oc (koR)^ 3 ^ 2 . Further depolarization inside one cell occurs, if Fara¬ 
day rotation by electrons around 7 min (Jones & O’Dell 1977) depolarizes the radiation 
within the cell 717 oc t/tf■ Here Tp is the Faraday depth in the cell, which is smaller 
than the Faraday depth for the whole source by a factor ( k^R ) _1 . The fractional linear 
polarization can then be approximated by 


7 T L ~ 


OL T 1 
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For sources above the IC limit the magnetic held strength must be very low and ra¬ 
diation at GHz-frequencies is emitted by electrons with large 7 ra( j. The resulting ratio 
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7 rad/ 7 min in Eq.[I] becomes large, and the polarization drops below 1% if 7 min < 10 3 at 
Tb = 10 12 K. At the eqnipartition temperature we still have depolarization if 7 min < 80. 
The fact that most IDV sources are variable in polarized flux with a mean 717 > 1% 
requires a cut-off in the electron population close to 7 , which dominates the radiation 
from r = 1 surfaces. At Tjj, > 10 13 K this requires hne-tuning of 7 min and provides strong 
constrains for acceleration mechanisms. Any substantial population of cold electrons 
with 7 min fa 1 is excluded in polarized sources with Tjj, > 10 10 K. This problem does not 
arise in pair plasma jets, because Faraday depolarization does not occur there. 


3 The alignment problem of thin shock propagation 


The thickness Az of a layer of post shock gas behind a relativistic shock, which travels 
with a shock-Lorenz-factor T 5 > 2 for a distance z is 


Az 



r l ~4 

r l -1 


(2) 


assuming that the gas leaves the shock at the sound speed and has an ultra-relativistic 
equation of state. A gas element that passed through the shock at z = 0 is separated 
by a distance Az from the shock, when the shock has travelled a distance z in the jet 
frame. For oblique shocks (as suggested by Spada et ah, 1999) the velocity and the 
thickness of the post shock gas will be larger. The ratio of thickness to the square-root 
of the surface area is 


r ~ Az/(y/nz sin 0 ) 



/(a/7t sin- 0 ) 


( 3 ) 


where it is assumed that the shock travelled at a constant velocity from the tip of a 
conical jet with half opening angle -0. In the following we will assume that the shock 
travels only 1/10 of that distance. When the shock is viewed face on, the surface area 
is a factor 1/f 2 larger than inferred from the variability timescale. But this factor is 
subject to relativistic aberration 


r = V sin 9 + T V r (cos 9 — (3 ) 


( 4 ) 


where [3 is the velocity of the post-shock gas and 9 the angle between jet and the line 
of sight. The true observed flux variations can arise from orientation changes of shock 
and jet 


AT obs 

F in t 


h = 3D 2 §r- 2 A« 


( 5 ) 


or from inhomogeneities in particle density or magnetic held along the jet (Qian et al., 
1991), which are scanned by the shock and travel with the sound speed through the 
post-shock gas 

= f 2 = D 3 (r/0.667)- 2 . (6) 

Uint 
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Figure 1: The square of the observable reduced timescale r^ 2 (dashed) and the Doppler 
boost T> ‘ > (dotted) for the brightness temperature in case of intrinsic variations are 
shown together with the flux-boosting factors f\ (solid) and f -2 (dash-dotted) as defined 
in the text. The misalignment is the angle 6 between jet direction and line of sight. 
The assumed shock is travelling with T^ = 6 in a jet of r pre = 2.7. The implied increase 
of observed Tj \ due to orientation changes f\ is not substantially larger than T> 3 and 
the result of scanned perturbations f<i is only sufficient for explaining IDV for strong 
alignment of jet and observer. 


In Fig.l the Doppler boosting of the flux due to intrinsic variations for a spherical 
source V 3 is compared to the flux-boosting factors f\ and f -2 for a shock with = 6 
in a jet with if = 2.5° and r pre = 2.7 for the preshocked gas (resulting in a post shock 
gas with T = 12). Orientation changes of thin shocks require A 6 ~ 6 and consequently 
fi can not substantially increase the flux variations on short timescales compared to 
variability in a spherical source. Only if disturbances along the jet are highlighted by 
the passage of a shock can the flux variations be so rapid and strong that /2 exceeds 
10 6 which is necessary to bring a source with intrinsic Tj, < 10 12 K up to the observed 
10 18 K. Nonetheless this cannot provide an explanation for most IDV sources, because 
it requires a misalignment of the line of sight to the jet direction much smaller than 
l/(3r), where the plateau of T> 3 starts (see Fig.l). In the specific case shown in Fig.l 
the misalignment must be less than lCV 3 r _1 , which cannot be reconciled with 25% 
of core-dominated sources showing IDV. Furthermore scanned perturbations in the jet 
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lead to variable optical depth in the shocked gas and timelags, which are not observed 
in cross correlations of IDV at different radio frequencies. 


4 Refractive Interstellar Scintillation 

It is known from pulsar measurements that compact radio sources are subject to scin¬ 
tillation in the ISM of our galaxy. Extragalactic sources flicker (Heeschen 1984) at 
frequencies v « 1 GHz with timescales of several days to weeks and this has been in¬ 
terpreted as a result of strong refractive scattering in the extended ISM of the galactic 
disk (Blandford, Narayan & Romani, 1986). The transition from strong to weak scat¬ 
tering is expected at about 5 GHz depending on the path length in the disk towards 
the source. The maximum of the modulation index m is ~ 1 for a point source. A lJy 
source with an apparent Tj [ of 10 12 K has a size of 230 /xas at A6cm and scintillation 
will be quenched[], because the source size is much larger than the angular size of the 
Fresnel scale in the ISM. For a typical distance of 200 pc to the scattering medium, 
the expected modulation index is m < 1.7%, while a source with T\\ = 10 13 K will have 
m < 5% with a characteristic timescale of 0.5 days. Therefore any lJy source with 
apparent T& above few xl0 12 K is expected to scintillate with a timescale comparable 
to those observed in IDV sources, as has been demonstrated by Rickett et al. (1995) 
for 0917+624. 

The light curves of scintillating sources contain further information, which can be 
explored via structure functions (SF) (Simonetti, Cordes & Heeschen, 1985): 

SF(T) = ([I(t+r)-I(t)] 2 ) i . (7) 

The SF extracted from the data can be compared to theoretical expectations. For 
strong refractive and weak scintillation in an extended medium Coles et al. (1987) 
gave a closed form for SF’s, which includes three cut-offs for the Fourier transform of 
the local SF. These are due to the Fresnel scale, the visibility amplitude, which contains 
the size of the source, and an exponential cut-off in strong scintillation giving rise to 
the refractive scale. Whichever comes first determines the modulation index and the 
slope of the SF below the first maximum (see Fig.2). Following Blandford, Narayan & 
Romani (1986) we assume a Gaussian distribution of ionized matter, keeping in mind 
that for IDV a much smaller scale height H ~ lOOpc is required than the ISM scale 
height in the galactic disk. 

We concentrate on the modulation index derived from the plateau of the SF at 
large timelags and on the slope for shorter lags. It turns out that in weak scintillation 
the slope of the SF is either quadratic SF(t) oc t a , a = 2 for steep turbulent spectra 
f3 > 4 or a = (3 — 2 for (3 < 4. 

In strong and quenched scintillation the SF shows a broken power law below the 
maximum and this slope also depends on the power-law index (3 of the turbulent 
spectrum <3>(g) = of density fluctuations in the ISM. For the limiting case f3 — 4 

1 A gaussian brightness distribution is assume for extended sources in this paper. 
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model 

H [pc] 

P 

C 2 n [m P 3 ] 

Ose [+as' 

6*520 [mas' 

v [km/s] 

K 

70 

11/3 

10~ 3 

40 

0.4 

5 

SI 

70 

4.567 

10- 10 - 6 

40 

0.16 

14 

S2 

140 

4.567 

10- 1215 

40 

0.16 

30 


Table 1: Parameters for the structure function models. The result for A 6 cm and A20cm 
are shown in Fig.2 and Fig.3. The angular size dse and $S 2 o are the source size at A 6 cm 
and A20cm. 

discussed in Blandford, Narayan & Romani (1986), this slope is a = 1 in the quenched 
and strong case, while for a Kolmogorov-likc spectrum (3 = 11/3 the slope is a = 5/6 
for strong scattering and a = 2/3 for quenched scintillation. For quenched scintillation, 
the slope is generally a = (3 — 3. 

We have computed theoretical structure functions dominated by quenched scintil¬ 
lation and compared them to the observed data taken at A20cm and A 6 cm (Fig.2 and 
Fig.3). The parameters are the power-law index /?, the strength of density fluctuations 
Cjf, the scale height H and the velocity of the observer relative to the ISM cloud. From 
simultaneous fitting of the SF’s at 6 cm and 20 cm we can find plausible parameters 
and the size of source at these wavelengths. The model parameters are summarised in 
Table 1. It is clearly visible from Fig.2 and Fig.3 that a steep spectrum is preferred. 

In particular the slope a ~ 1.3 at A20cm indicates a steep spectrum (3 = 4.3, when 
only quenched scintillation is assumed. The turbulent spectrum has to be even steeper 
because the SF flattens before turning over into the a = (3 — 2 slope dictated by the 
Fresnel cut-off at smaller timelags. Together with the weak curvature near the plateau 
for large timelags, this points to the slope (3 ~ 4.6 of the models SI and S 2 . 

5 Conclusions 

Based on IDV observations of several Quasar and BL Lac radio cores with apparent 
brightness temperatures in the range of 10 16 -10 21 K we investigated the possibility of 
intrinsic variations due to the propagation of thin relativistic shocks. We find that no 
source with 7/ > 10 16 K can be explained by that model without calling for Doppler 
factors larger than 20 and strong alignment between the jet and the line of sight. 
The required alignment firmly rules out this hypothesis as an explanation for all IDV 
sources. 

Based on refractive interstellar scintillation IDV can be explained by turbulence in 
the ISM. The scale size of the ionized gas responsible for IDV is about 100 pc in the 
case of 0917+624, and scintillation is quenched by the source size, which is larger than 
the refractive or Fresnel scale in the ISM. A fit to the structure functions at A20cm 
and A 6 cm indicates a steep power law (/3 ~ 4.6) for <L(g) fluctuations. This corre¬ 
sponds to an energy spectrum ~ q~ 2 - 6 that is much steeper than both a Kolmogorov 
spectrum and the q~ 2 spectrum for compressible turbulence, but a bit shallower than 
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0917+624 at 20cm 



Figure 2: Structure function of total flux variations of 0917+624 at 20 cm in 1989. 
The data have recently been analysed by Qian et al. (2001). Over-plotted on the SF 
derived from the measurement are models for quenched scintillation according to Table 
1. The dash-dotted line shows the model with a Kolmogorov spectrum (K). The steep 
spectra models (SI - dashed) and (S2 - solid) give better fits for the 20cm data. 


~ q~ 3 expected for two-dimensional turbulence. The length scale probed in the ISM 
by these measurements are between 5 10 8 m and 2 10 9 m at the peak of the structure 
functions. The slope of the turbulent spectrum is derived from structure functions at 
small timelags and the corresponding spectrum extends an order of magnitude down 
to smaller spatial scales. Compressible turbulence is not unexpected at these scales, if 
turbulence is driven by shocks from supernovae or by stellar winds. 

Scintillation in 0917+624 is quenched by the source size, which is one parameter of 
the theoretical fits to the SFs. The required sizes are ~ 40/ias and ~ 0.4 mas at A6cm 
and A20cm respectively. At A6cm, the flux of 0917+624 (redshift z = 1.446) is 1.5Jy. 
Combined with the derived source size this implies 77 = 10 14 K. Again Doppler factors 
of about 100 are needed to avoid the IC catastrophe or Doppler factors of 1000 to arrive 
at equipartition temperature. The sizes derived for 0917+624 from structure function 
models might be changed, if a degeneracy in the model parameters exists. The most 
plausible direction is an even closer ISM screen with a higher level of turbulence. In 
this case the angular Fresnel scale gets larger and larger source sizes are allowed. 
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0917 + 624 at 6cm 



Figure 3: Structure function of total flux variations of 0917+624 at 6cm in 1997. The 
data come from one of the best sampled light curve at this wavelength. Over-plotted on 
the SF derived from the measurement are the same models for quenched scintillation 
like in Fig. 2 applied to A6cm. 
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